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Abstract. We use a novel statistical tool, the mark correlation functions (MCFs), to study clustering of galaxy- 
size halos as a function of their properties and environment in a high-resolution numerical simulation of the 
ACDM cosmology. We applied MCFs using several types of continuous and discrete marks: maximum circular 
velocity of halos, merger mark indicating whether halos experienced or not a major merger in their evolution 
history (the marks for halo with mergers are further split according to the epoch of the last major merger), and 
a stripping mark indicating whether the halo underwent a tidal stripping (i.e., mass loss). We find that halos 
which experienced a relatively early (z > 1) major merger or mass loss (due to tidal stripping) in their evolution 
histories are over-abundant in halo pairs with separations < 3/i _1 Mpc. This result can be interpreted as spatial 
segregation of halos with different merger histories, qualitatively similar to the morphological segregation in the 
observed galaxy distribution. In addition, we find that at z — the mean circular velocity of halos in pairs of halos 
with separations < 10/i _1 Mpc is larger than the mean circular velocity U c irc of the parent halo sample. This mean 
circular velocity enhancement increases steadily during the evolution of halos from z = 3 to z = 0, and indicates 
that the luminosity dependence of galaxy clustering may be due to the mass segregation of galactic dark matter 
halos. The analysis presented in this paper demonstrate that MCFs provide powerful, yet algorithmically simple, 
quantitative measures of segregation in the spatial distribution of objects with respect to their various properties 
(marks). This should make MCFs very useful for analysis of spatial clustering and segregation in current and 
future large redshift surveys. 
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1. Introduction 

The advent of large wide-field redshift su rveys of galax- 
ies, such as the Two-Degree Field (2dF, Colless et al. 



2001) and the Sloan Digital Sky Survey (SDSS, York ct 



luminosity (se e, e.g., Hamilton 1988 ; Bcnoist et al. 1996 ; 
Norberg 2001 ), and colors (e.g., Zehavi, 2002 ) of galaxies. 
Greatly enhanced clustering of super-luminous IR-galaxies 



al., 20Q€), will allow detailed sLudies uf clusLcihig uf galax- 



( Bouchet et al., 1993 ) and morphological ( Dressier, 1980 



ies as a function uf llieh eiivhuuuieuL and iiiLeiiial piup- 

erties. Indeed, hierarchical growth of structure via gravi- 
tational instability is thought to play a dominant role in 
shaping both the large-scale galaxy clustering and internal 
properties of galaxies such as luminosity, colors, and mor- 
phology. This close connection implies that studies of the 
spatial distribution of galaxies as a function of their inter- 
nal properties and environment should provide us valuable 
insights into the process of galaxy formation. Previous 
observational studies and the first results from the 2dF 



et al., 2002|) and color (e.g., [Butcher fc Oemler, 1978 



Postman fc Gcllcr, 1984| ; |Whitmorc ct al., 1993^ [Bivianc 



Margoniner et al., 2001 ) segregation in clusters of galaxies 
indicate dependence of clustering on the merging history 
and large-scale environment. 

In this paper we use a novel statistical tool, the mark 
correlation functions (hereafter MCFs), to study clus- 
tering of galactic halos as a function of their proper- 
ties and environment in a high-resolution numerical sim- 
ulation of the ACDM cosmology. Mark correlation func- 



and SQ SS have oh own that cluotoring otrongth depondo o n 



tions (Stoyan, 1984; Stoyan & Stoyan, 1994 ) have been 
introduced into astrophysics only recently ( Beisbart fc 



morph ology (e.g., | Hormit ct al., 1096 | ; | Cuzao ct al., 1997| ) 



Kerscher, 2000 ), although so me aspects of m arked point 
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processes were discussed by Peebles (1980 ). The mark 
statistics can be used to quantify the differences in the spa- 
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tial distributions of various galaxy samples (similarly to the four year COBE DMR observations (Bunn & White 



the usi al two-point correlation function) and, at the same 



time, to study the interplay between the spatial clustering 
and the distribution of galaxy properties (marks). In this 
respect, the MCFs are the natural extension of the spatial 
correlation functions for studies where it is advantageous 
to consider mark and spatial distributions simultaneously. 
Variants of the MCFs can be used to study both continu- 
ous (e.g., luminosity or angular momentum) and discrete 
distributions (e.g., color classes or morphological types) of 
galactic properties. This makes them valuable for quanti- 
tative studies of luminosity and morphological segregation 
of galaxies as well as dependence of spatial distribution 
of galaxies on various events in their evolutionary history 
(e.g., time since the last major merger), which can be used 
as marks. Indeed, the mark correlation statistics quickly 
proved to be a very useful tool for identification of physi- 
cal processes that shape the observed spatial distribution 
of galaxies ( Szapudi et al. 200C , see Beisbart et al. 2002 
for a recent review) . 

The first step towards the use of clustering as a probe 
of processes shaping properties of galaxies is a good the- 
oretical understanding of how these processes affect spa- 
tial distribution of galaxies. During the last several years, 
thanks to continuously improving spatial and mass reso- 
lutions of numerical simulations and development of semi- 
analytic models of galaxy formation, there was a signifi- 
cant progress in the theoretical understanding of galaxy 



clustering evolution and bias (e.g., 


Bagla, 1998|; 


Jing, 


1998; 


<[auffmann et al., 1999; Katz et al., 1999|; 


Colin 


et al., 


999; 


Kravtsov & Klypin, 1999; 


Pearce et al., 1999; 


Schmalzing et al., 1999 


). The complete information about 



the internal properties and evolution of galactic halos, usu- 
ally available in theoretical analysis, allows one to study 
in detail the interplay between different evolutionary pro- 
cesses and spatial distribution of objects. In the present 
paper we use mark correlation functions to study cluster- 
ing of galaxy-size dark matter halos and its dependence 
on the halo properties (e.g., circular velocity) and evo- 
lution history in a high-resolution simulation of the cur- 
rently favored flat Cold Dark Matter (CDM) model with 
cosmological constant (see § ||). 

The paper is organized as follows. We discuss briefly 
the numerical simulation in Sect. |[ In Sect. || we introduce 
and explain the properties of mark correlation functions. 
In Sect. |^ we present analysis of the spatial distribution 
of dark matter halos using mark correlation functions. In 
Sect. |5| we discuss results and their implications and sum- 
marize our conclusions. 

2. Numerical simulations 

We used the Adaptive Refinement Tree (ART) code 



(Kravtsov et al., 1997) to simulate the evolution of colli- 
sionlcss DM in the currently favored ACDM model (il m = 
1 - fi A = 0.3; H = h x 100km s" 1 Mpc" 1 = 70km s" 1 
Mpc -1 ; erg = 1.0). The age of the Universe in this cosmol- 
ogy is 13.5 Gyrs and normalization is in accordance with 



1997) as well as the observed abundance of galaxy clusters 
(e.g., [Pierpaoli et al., 200l| ; [Ikebe et al., 2002j ). 



The numerical simulation of the ACDM model fol- 
lowed the evolution of 256 3 « 1.67 x 10 7 particles in 
a periodic 60/i _1 Mpc box. The particle mass is thus 
w 1.1 x 10 9 /i _1 M Q . The ART code reaches high force res- 
olution by refining all high-density regions with an auto- 
mated refinement algorithm. The refinements are recur- 
sive: the refined regions can also be refined, each subse- 
quent refinement having half of the previous level's cell 
size. This creates an hierarchy of refinement meshes of 
different resolution covering regions of interest. The crite- 
rion for refinement is local overdensity of particles: in the 
simulation presented in this paper the code refined an in- 
dividual cell only if the density of particles (smoothed with 
the cloud-in-cell scheme) was higher than n t h — 5 parti- 
cles. Therefore, all regions with overdensity higher than 
S = n t h 2 31 /n, where n is the average number density of 
particles in the cube, were refined to the refinement level 
L. For the simulation presented here, n is 1/8. The peak 
formal dynamic range reached by the code on the high- 
est refinement level is 32, 768, which corresponds to the 
smallest grid cell of 1.83/i _1 kpc; the actual force resolu- 
tion is approximately a factor of two lower. The simulation 
that we analyze here has been used in |Cohn et al. (1999 ), 
and we refer the reader to this paper for further numerical 
details. 

Identification of DM halos in the very high-density en- 
vironments (e.g., inside groups and clusters) is a challeng- 
ing problem. The goal of this study is to investigate spatial 
correlations of halo populations as closely related to the 
observed galaxy population as possible. This requires iden- 
tification of both isolated halos and satellite halos orbit- 
ing within the virial radii of larger systems. The problems 
associated with halo identification within high-density re- 



gions are discussed in Klypin et al. (1999). In this study we 
use a halo finding algorithm called Bound Density Maxima 
(BDM). The main idea of the BDM algorithm is to find 
positions of local maxima in the density field smoothed 
at a certain scale and to apply physically motivated cri- 
teria to test whether the identified site corresponds to a 
gravitationally bound halof|. It is based on the idea s of the 
DENMAX halo finder flBertschinger fc Gelb, 199l] ), in the 



sense that the BDM makes sure that the density peaks are 
gravitationally bound and estimates parameters of the ha- 
los after removing unbound particles. The algorithm iden- 
tifies both isolated halos and subhalos located in the virial 
regions of more massive halos. The distribution of halos 
identified in this way can be compared to the distribution 
of galaxies directly, because the halo and galaxy catalogs 
include both isolated systems and objects within clusters 
and groups. 

Even with algorithms tailored for identification of sub- 
halos, additional problems exist. Interacting halos ex- 



1 The detailed d escript ion of the algor ithm is given in 
( |Klypin et al., 1999| ) and ( |Coh'n et al., 1999| ). 
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change and loose mass; the total mass of a halo depends 
on its radius, which is difficult to define in a dense envi- 
ronment within virialized regions. We alleviate the latter 
problem by using the maximum circular velocity instead 
of the mass. In practice, maximum circular velocity is a 
rather stable quantity which changes little even when ha- 
los looses most of its mass and can serve therefore as a 
useful mass-related "tag" of a halo. Numerically, the maxi- 
mum "circular velocity" (y/GM/R), v c i TC can be measured 
more accurately then mass. In addition, the maximum cir- 
cular velocity can be more readily compared to observa- 
tions than, for example, virial mass or mass within the 
tidal radius of the halo. 



3. Mak correlation functions 



In studying galaxy clustering with the mark correlation 
functions, we view galaxies as discrete points in space with 
marks describing their intrinsic physical properties. Thus, 
we consider a point set {x.i}fL 1 and attach a mark m, 
to each point Xj £ K 3 resulting in the marked point set 



{(x„m,)}f =1 QStoyan, 1984) ; |5toyan fc Stoyan, 1994| ). The 
marks, in turn, can be either continuous or discrete. In 
the following, we use the circular velocity as a continu- 
ous mark and merging/stripping events of halos as dis- 
crete marks. In a subsequent paper, we will apply MCFs 
to study various other marks, such as the spin parame- 
ter (continuous scalar mark) and the angular momentum 
(vector mark). 

Let g be the mean number density of the points in 
space and g M (m)dm the probability that the value of a 
mark on a point lies within the interval [m, m + dm]. The 
mean mark is then m = f dm g M (m)m, the mark variance 
is cr| f = J dm g M (m)(m — m) 2 . We assume that the joint 
probability g SM (x, m) of finding a point at position x with 
mark M, splits into a space-independent mark probabil- 
ity and the mean density: g M (m) x g. The spatial-mark 
product-density 



£*2 M (( x i> TO i): ( x 2,m 2 )) dVxdm! dV 2 dm 2 



is the joint probability of finding a point at Xi with the 
mark m\ and another point at X2 with the mark m 2 . We 
obtain the spatial product density g>2(xi,X2) and the two- 
point correlation function £(r) by marginalizing over the 
marks: 



g 2 (1 + eW) = ea(xi,x a ) 



dmi I dm 2 g 2 M ((xi, mi), (x 2 , ma)), (2) 



where £(r) is the spatial two-point correlation function 
which depends only on the separation r = |xi — x 2 | of the 
points for a homogeneous and isotropic point set. 
We define the conditional mark density: 



A4 2 (mi,m 2 |xi,x 2 ) 



gf M ((xi,mi),(x 2 ,m 2 )) 



2 (xi,x 2 ) 



for g 2 (xi,x 2 ) ^ 0, 
else . 



(3) 



For a stationary and isotropic point distribution 
■A4 2 (mi, m 2 |r)dmidm2 is the probability of finding the 
marks mi and m 2 at two points located at xi and x 2 , un- 
der the condition that they are separated by r = |xi — x 2 |. 
Now the full mark product-density can be written as 

£>f M ((xi,mi), (x 2 ,m 2 )) = A4 2 (mi,m 2 |xi,x 2 ) £> 2 (xi,x 2 ). 

(4) 

If there is no mark segregation in space, M. 2 (m\ 1 m 2 \r) is 
independent of r, and Ai 2 (mi,m 2 \r) — g M (mi) g M (m 2 ) . 

Starting from these definitions, especially using the 
conditional mark density M 2 (mi,m 2 \r), one may con- 
struct several mark correlation functions sensitive to dif- 



ferent aspects of mark segregation (Bcisbart & Kerscher, 



2000). The basic idea is to consider weighted conditional 



correlation functions describing the probability of find- 
ing two points at a separation r. For a positively defined 
weighting function /(mi, 7712) > we define the average 
over pairs with separation r: 



(/)p ( r ) = / dmi J dm 2 f(mi,m 2 ) M 2 (m u m 2 \r). 



(5) 



(/)p ( r ) i s the expectation value of the weighting function 
/ (depending only on the marks) , under the condition that 
we find a point pair with separation r. For a suitable de- 
fined integration measure Eq. (||) is also applicable to dis- 
crete marks. The definition (|5|) is very flexible, and allows 
us to investigate the correlations of both continuous and 
discrete marks. 

In the following analysis, we calculated the mark cor- 
relation functions taking into account the periodicity of 
the simulation box. However, we obtain virtually identical 
results using the estimator without boundary corrections 
(see Appendix A of Beisbart fc Kerscher (2000 ) for de- 
tails). 

3.1. Correlations of scalar marks 



(1) For scalar marks the following mark correlation functions 



have proven to be useful (Stoyan & Stoyan 1994, Beisbart 
|fc Kerscher 2000| , |Schlathcr 2001] ): The simplest weight to 
be used is the mean mark: 



k m (r) 



(mi + m 2 ) p (r) 
2 m 



(6) 



It quantifies the deviation of the mean mark on pairs with 
separation r from the overall mean mark m. For example, 
km (f) > 1 indicates mark segregation for point pairs with 
a separation r, specifically their mean mark is larger than 
the overall mark average. 

Higher moments of marks like the mark fluctuations 



var(r) = ( (mi — (mi) P (r)) ) (r) 



(7) 



or the mark covariance ( Cressie, 1991 ) 

cov(r) = ((mi - (mi)p (r)) (m 2 - (m 2 } P (r))) p (r) 

= (mim 2 } P (r) - (mi) p (r) (m 2 ) p (r), (8) 
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may be used to quantify mark segregation. For example, a 
positive ccw(r) indicates that points with separation r tend 
to have similar marks, whereas a negative cov(r) indicates 
different marks. 

3.2. Correlations of discrete marks 

For discrete labels only combinations of indicator- 
functions are possible, and the integral degenerates into 
a sum over the labels. Supposing the marks of our ob- 
jects belong to classes labeled with the conditional 
cross-correlation functions are given by 

Cij(r) = (S mi iS m2 j + (1 - Sij)5 m2i 5 mi j) p (r), (9) 

with the Kroncckcr 5„ u j = 1 for mi = i and zero otherwise 
flBtoyan fc Stoyan, 1994|) . By construction 

EE^w = L ( 10 ) 

i 3 

Mark segregation is indicated by Cy- ^ %QiQj/ Q 2 for i ^ j 
and Cn 7^ gf j ' g 2 , where gi denotes the number density of 
points with label i. The C« are cross-correlation functions 
under the condition that two points are separated by a 
distance of r. 
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Fig. 1. Evolution of the cumulative number of halos with 
a circular velocity larger than v c irc- The solid, dot-dashed, 
and dashed curves correspond to z = 0, z = 1, and z = 3, 
respectively. 



4. Results 

For any study one needs to have a complete halo sample 
that is not affected by the numerical details of the halo 
finding procedure. We have tested the completeness of the 
halo samples using different parameters for the halo finder. 
For the given force and mass resolution the halo samples 
do not depend on the numerical parameters of the halo 
finder for halos with v c i rc ^100 km/s ( Gottlober et al 
1999| ). In Fig. we show the cumulative number of halos 



with a circular velocity larger than a value w c i rc , for red- 
shifts z = 0, z = 1, and z = 3, respectively. Assuming 
a minimum circular velocity of 100 km/s the samples are 
complete at z < 1 but we are missing a small fraction of 
halos with w C j rc < 130 km/s at z = 3. 

4.1. Merging of halos 

According to hierarchical structure formation halos 
formed early and grow during evolution due to accretion of 
matter and merging with other halos. In particular merger 
events are important because they are expected to lead to 
dramatic changes in the structure of dark matter halos and 
the galaxies they harbor. In-falling objects may damage 
or even destroy a stellar disk. The inflow of material may 
also serve as a source of fresh gas and therefore increase 
the star formation rate. At the same time, collisions be- 
tween halos may result in shock heating of the gas, which 
would tend to delay or prevent star formation for some 
period of time. 



than 25% during a time interval of about 0.5 Gyrs (ap- 
proximate interval between simulation outputs). In the 
above paper we showed that for redshifts z < 2 the merger 
rate can be fitted by a simple power law (1 + z) 3 . This 
merger rate evolution is in very good general agreement 



Following Gottlober et al. (2001 ) we identify major 
mergers as events when the mass of a halo grows by more 



with observations (e.g., |Le Fevre et al. (2000 ) measured 
a merger rate varying with redshift as oc (1 + z) 3 2±0 - 6 ). 
In addition, we found that evolution of the merger rate 
depends on the environment of the halo: halos that end 
up in clusters and groups by z = have a steeper evo- 
lution of merger rate and a higher rate of major mergers 
at early epochs compared to isolated "field" halos. This is 
because clusters and groups form in the regions that are 
overdense on large scales in which halos form and undergo 
the phase of active merging earlier than the overall field 
halo population. 

For the z = sample of halos with i; c ; rc > 100 km/s, 
about 32 % of halos had one major merger in the past and 
additional 19 % of halos had two or more major mergers. 
Now let us consider the distribution of epochs of the last 
major merger (relevant, for instance, for estimating a frac- 
tion of halos that could host old disks such as that of the 
Galaxy). We found that 55 % of the halos located in clus- 
ters at z — underwent a major merger after redshift 
z = 4, but that corresponding fraction for the isolated 
halos is 43 %. In contrast, the fraction of isolated halos 
which underwent a major merger at a redshift z < 1 is 
somewhat higher (19 %) than the corresponding fraction 
of halos in clusters (14 %); for z < 0.4 (i.e., within the 
last « 5 Gyrs) these fractions are 8 % and 3.5 %, respec- 
tively. This reflects the fact that due to the high internal 
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velocity dispersion of halos in clusters mergers are almost 
impossible. Since the merger rate of group halos is high 
compared to the overall merger rate of halos at all an- 
alyzed redshifts, the fraction of present-day group halos 
which merged after a fixed redshift is always higher than 
the fraction of isolated or cluster halos which merged after 
the same redshift. 

Finally, let us consider one more effect. Due to the tidal 
interactions halos in dense regions (i.e., virialized regions 
of groups and clusters) tend to loose mass via tidal strip- 
ping. In order to take this effect into account, we follow 
the mass evolution of all halos and identify halos that lost 
more than 30% of their maximum (over their evolution) 
mass from z = 1 to the present epoch. One would expect 
that galaxies hosted by such halos also lost their supply 
of fresh gas so that no star formation was possible in the 
recent past. 

4.2. Spatial distribution and evolution history 

In the preceding section we considered the overall fractions 
of halos in different environments and with certain merger 
history classes. This relatively straightforward analysis re- 
veals existence of some environmental dependency of halo 
evolution histories. The goal of this section is to carry 
out a more quantitative analysis of how spatial clustering 
of different halo subsamples depends on evolution histo- 
ries of their halos. As discussed in the previous subsection, 
major mergers (and tidal stripping) can be expected to re- 
sult in dramatic changes in the properties of galaxies (i.e., 
morphology and color). One can expect, therefore, that 
the spatial distribution of halos that experienced a recent 
merger or stripping event is different from the distribu- 



tion of the overall halo population. For example, Knebe & 



Muller (1999) found that massive halos undergoing merg- 
ers at present exhibit a much stronger bias with respect 
to the dark matter than relaxed halos do. 

Figure [2] shows the two-point correlation function of all 
halos (a) with w c ; rc > lOOkm/s (solid line) compared with 
that of the subsamples of halos with different evolution 
histories. We divided the sample of all halos into four sub- 
samples: halos which never (n) underwent a major merger 
in the past, halos which underwent a major merger before 
(e) and after (1) redshift z = 1, and halos that lost mass 
since z — 1 (s). Note, that the stripped halos constitute 
a separate sample, however every stripped halo belongs 
also per definition to one of the other subsample. In par- 
ticular, a substantial part of "stripped" halos in clusters 
underwent a major merger before redshift z = 1, i.e. they 
belong to the sample (e) of halos. They lost most of their 
mass later on due to interactions. 

Figure |^ shows clearly that the subsample of halos 
which underwent the last major merger before redshift 
z = 1 is more clustered than the sample of all halos. This 
is not surprising since most of such halos formed early 
in the regions of large-scale overdensity and ended up in 
groups and clusters by the present epoch. The stripped ha- 




r [h 1 Mpc] 



Fig. 2. The two-point correlation function of all (a, solid 
line) and subsamples of halos with different mass evolution 
histories (e: last major merger at z > 1, long-dashed line; 
1: last merger at z < 1, short-dashed line; n: never merged, 
dot-dashed line; and s: stripped halos, dotted line. 



los are even more biased with respect to the overall halo 
population. This is also due to the fact that halos that 
loose mass via stripping are located in the high-density 
regions where tidal forces are strong. 

Let us now consider the spatial distribution of dif- 
ferent halos using the mark correlation functions intro- 
duced in Sect. ||. To this end, we first split the total halo 
sample (i.e., sample of all halos with circular velocities 
^circ > lOOkm/s) into two subsamples consisting of ha- 
los which experienced a major merger (sample m) and 
halos which never experienced a major merger (sample 
n), respectively. Figure |^ shows the conditional cross- 
correlation functions (eq. [)]) of these samples at z = Q. 
Positive C m , m on scales below 3/i _1 Mpc indicates that 
halos that experienced a major merger in their formation 
history are relatively overabundant in close pairs of halos, 
while halos without a major merger are underabundant. 
No significant cross-correlation between m and n exists. 
Qualitatively, this feature is independent from a lower cut 
in the circular velocity i> c i rc , but the amplitude of the ef- 
fect is reduced if we consider only the more massive halos 
with v c i rc > 150km/s. To investigate the significance of 
these deviations we performed a non-parametric Monte 



Carlo test ( 


Besag & Digglc, 1977 


) , similar to the 


in Kcrscher ct al. (2001 


) and 


Kerscher (1998). 



hypothesis is "no mark correlation" . We simulate this null 
hypotheses by keeping the positions of the halos fixed and 
randomly re-assigning the marks (in this case the class 
assignments) to the halos. We repeat this to generate 99 
realizations of this null hypothesis. The shaded areas in 
Fig. [| are the one-c regions numerically determined from 
these samples. To quantify the significance we have to de- 
fine a distance measure. Using M = 4 equidistant radii ri 
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Fig. 3. Conditional cross-correlation function of halos 
with (m) and without (n) major merger in their evolu- 
tion history. The shaded area is obtained by randomizing 
the assignment among the halos (see text for details) . The 
subsamples were selected from the total sample of halos 
with t> c j rc > lOOkm/s. 

in the range from 0.8/i _1 Mpc to 2.8/i _1 Mpc we define the 
"distance" of the k-th sample to the expected value for 
no mark correlation: 



d k = 



1 

M 



M 



0^ 



(11) 



with the overall number density g and the number density 
g m of merged halos. Similarly we determine the distance 
d hal ° of the original halo sample to the null hypothesis. 
Then we sort all the d's and determine the position of 
d hal °. In this case d hal ° is the fifth largest distance, and 
we conclude that the original halo sample is incompatible 
with the null hypothesis "no mark correlation" at a signif- 



icance level of 95% (see the comments by Marriott 1978 
concerning the significance level). 

Figure || shows the conditional cross-correlations of 
the three subsamples which we considered above using the 
two-point correlation function (Fig. ^: the halos which 
never underwent a major merger in the past, sample n; 
early major merger at z > 1, sample e; and late major 
merger at < z < 1, sample 1). 

As before, positive C Cj0 at small separations indicates 
that for objects at distances less than about 2 h~ 1 Mpc the 
halos with a major merger in their early formation history 
are relative over-abundant, at the expense of halos with- 
out a major merger, as deduced from the lowered C njI1 . 
This signal is most prominent on scales below 3/i _1 Mpc 
but it extends out to 10/i _1 Mpc in agreement with the en- 
hanced two-point correlation function of that subsample 
(Fig. H). We interpret this as indication of a high num- 
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Fig. 4. Conditional cross-correlations of halo subsam- 
ples defined using the epoch or absence of the last major 
merger (e: the last major merger at z > 1, 1: at z < 1, n: 
never occured) . The shaded area is obtained by randomiz- 
ing the assignment among the halos (see text for details). 
The subsamples were selected from the total sample of 
halos with v C i rc > lOOkm/s. 



ber of early merged halos in clusters. The signal has high 
significance and it is not influenced by uncertainties in 
the normalization of the correlation function that may be 
caused by selection effects. Interestingly, the halos with a 
late major merger show no excess correlations but rather a 
lowered abundance on small scales, also manifested as the 
lower correlation function amplitude of that subsample. 

The lower panel of Fig. |] shows the cross-correlation of 
halos from different evolution classes. The over-abundance 
of pairs of never and early merged halos reflects the contin- 
uous accretion process onto high density regions. Inf ailing 
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isolated halos from less dense regions accrete onto higher 
density regions with high velocity dispersions and, thus, 
low probability of merging. Therefore, type-n halos can 
survive in the high-density regions relatively long which 
explains the increasing of C nje towards small scales. The 
opposite is true for never n and late I merged halos. Type- 
1 halos located predominantly in groups where mergers 
are more likely due to the lower velocity dispersions. The 
probability of accreting type-n halo (located close to an 
type-1 halo) to experience a merger in high-density regions 
is therefore high. Many such halos will thus disappear (will 
become 1-halos) resulting in suppression of C ni i amplitude 
at small separations. 

Note, that these features are qualitatively indepen- 
dent from a lower cut in the circular velocity v c \ TC , but 
the amplitude and the spatial range is reduced if we 
consider more massive halos with w c i rc > 120km/s or 
w c irc > 150km/s. This is due to the higher number of 
mergers within the low circular velocity halos. 

Let us now consider halos which lost a substantial part 
of their mass due to tidal interactions. Figure || shows the 
conditional cross-correlation function of the halos using 
two classes: no stripping (sample: ns) and stripped (sam- 
ple: s). Positive C SjS at r < 5/i _1 Mpc indicates that the 
number of stripped pairs is strongly enhanced at these 
separations, whereas the number of non-stripped pairs is 
reduced. This result is in accordance with the strongly en- 
hanced correlation function shown in Fig. ||. In fact, we 
expect to find stripped halos only in the environment of 
clusters. The results for samples with a higher cut in the 
circular velocity v c i rc > 150km/s are very similar. 



4.3. Spatial distribution and circular velocity 

As mentioned in Sec. ^, the mass and the maximum circu- 
lar velocity of halos are tightly related. At the same time, 
the circular velocity can be determined more reliably in 
simulation as well as in observations, either through di- 
rect measurement using emission line width or rotation 
curve or via galaxy luminosity using the Tully-Fisher and 
the Faber-Jackson relations. Therefore, it is interesting to 
explore galaxy mark correlations with the maximum cir- 
cular velocity as mark. This would mimic to some degree 
luminosity segregation effects in observed galaxy samples. 
Figure || shows the mark correlation functions of halos at 
z = with the circular velocity as mark. There is a strong 
signal in k m (r) at small separations but the signal is signif- 
icant even out to 10/i _1 Mpc . This indicates that the mean 
circular velocity of pairs of halos with separations below 
~ 10ft, _1 Mpc is larger than the overall mean circular ve- 
locity ?J c i rc of the parent halo sample. The negative signal 
of cov(r) is confined to small separations (< 2/i _1 Mpc). 
This is because pairs at small separations are more fre- 
quently built from one halo with circular velocity larger 
than TJcirc and the other halo with circular velocity smaller 
than zT C i lc . Hence, this signal is dominated by pairs of 
small-mass subhalos and massive parent halos. The mark 
correlations results for both samples with w circ > lOOkm/s 
and u c irc > 150km/s are shown. The mean mark correla- 
tion function, k m (r), exhibits a slightly stronger signal for 
the sample with a lower cut in circular velocity, but the 
signal cov(r) is weaker for t> circ > lOOkm/s. The latter is 
due to the considerably larger number of isolated small- 
mass halos; i.e. in addition to pairs between large- and 
small-mass halos, for the u c irc > lOOkm/s sample there 
are many more small-small mass pairs. 

Figure ^ shows the evolution of the conditional covari- 
ance cov(r) of the circular velocity with redshift. The con- 
ditional covariance is negative at low redshifts out to scales 
of 2/i _1 Mpc as discussed above. At high redshifts (z = 2 
and z = 3), significant positive amplitude of cov(r) indi- 
cates that pairs with similar circular velocities are over- 
abundant. At these redshifts the signal is significant out 
to the scale of ~ 5/i _1 Mpc due to the large number of 
smaller- mass progenitors of the present day halos. 

5. Discussion and conclusions 

In the previous section we used a novel statistical tool, the 
mark correlation functions, to study clustering of galactic 
halos as a function of their properties and environment in 
a high-resolution numerical simulation of the ACDM cos- 
mology. We applied MCFs using several types of continu- 
ous and discrete marks: maximum circular velocity, u c i rc , 
of halos (continuous) , merger mark indicating whether ha- 
los experienced (m) or not (n) a major merger in their 
evolution history, a stripping mark (s) indicating whether 
the halo underwent a tidal stripping (i.e., mass loss) dur- 
ing its evolution (discrete marks). The halos which under- 
went major merger (m) are further classified by the epoch 
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Fig. 6. The mark correlation functions k m (r) and cov(r) 
(see Eqs. [| and ||) for halos with the circular velocity 
w c irc as a scalar mark. Two lower cut-offs have been used: 
Wcirc > lOOkm/s (solid line) and u c i rc > 150km/s (dot- 
ted line). The shaded area is obtained by randomizing the 
mark among the halos. 

of the merger: late (1; z < 1) and early (e; z > 1) mergers. 
Our main results are as follows. 

The two-point correlation amplitude is different for the 
halo subsamples with different marks. The halos that ex- 
perienced an early major merger or mass loss (e and s) are 
clustered considerably more strongly than the overall halo 
population, while halos with late or no mergers have cor- 
relation function amplitude below that of the overall halo 
sample. This result indicates that halo clustering depends 
sensitively on the details of their evolution history. If exis- 
tence of a major merger during halo evolution is related to 
the morphology of galaxies that halos host, the above re- 
sult indicates that early type galaxies and galaxies in clus- 
ters and groups (hosted by halos that undergo tidal strip- 
ping) should be clustered more strongly than the late type 
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Fig. 7. The redshift evolution of the mark correlation 
function cov(r) for halos with the circular velocity w C i rc 
as a scalar mark and a lower cut-off of i> c j rc > lOOkm/s: 
z = (solid line), z = 1 (dotted line), z — 2 (dashed line), 
z = 3 (dashed dotted line). The shaded area is obtained 
by randomizing the mark among the halos. 



galaxies and the overall galaxy population. Qualitatively, 
such trend exists in the observed galaxy samples (e.g., 
Hermit et al., 1996| ; |Guzzo et al., 1997| ; |Zchavi, 2002j) 



im- 
plying that the morphology-dependent clustering may be 
largely due to the overall merger history of the galactic 
halos. 

Using maximum circular velocity of halos as a con- 
tinuous mark, we found that at z = the mean circular 
velocity of pairs of halos with separations < 10/i _1 Mpc 
is larger than the overall mean circular velocity u C irc of 
the parent halo sample (manifested as significant enhance- 
ment of the mean mark at these separations; see Fig ^|). 
Moreover, the negative mark covariance (Eq. ||) at small 
separations shows an enhanced abundance of pairs with 
halo circular velocities above and below the average cir- 
cular velocity. This mean circular velocity enhancement 
increases steadily during the evolution of halos from z = 3 
to z = 0. The mark covariance, cov(r), has a more com- 
plicated behavior: it is negative at present, disappears at 
redshift z ~ 1 and becomes positive at higher redshifts 
due to the larger number of low circular velocity halos 
(the circular velocity function of halos steepens at high 
redshifts, see Fig. |l|). Although the relation is not direct, 
the maximum circular velocity of halos should correlate 
well with the maximum circular velocity or velocity dis- 
persion of the galaxies they host. The enhanced mean cir- 
cular velocity in small-separation pairs should therefore 
correspond to the luminosity segregation or luminosity- 
dependent clustering in the observed galaxies. The lumi- 
nosity dependence of galaxy cl ustering was re cently con- 
vincingly dete cted in both 2dF flNorberg, 200l|) and SDSS 
( Zchavi, 2002 ) galaxy surveys. 
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The mark correlation analysis indicate that galaxy- 
size halos (w c irc > lOOkm/s) which experienced a major 
merger in their evolution history are over-abundant in 
pairs with separations < 3 /i Mpc with respect to the 
overall halo population, while halos which never experi- 
enced a major merger are under-abundant at these sep- 
arations. We find no significant cross-correlation between 
these two halos classes. The overabundance of merger ha- 
los is due largely to the halos which experienced a major 
merger relatively early (z > 1); halos with late (z < 1) 
major merger are not over-abundant (this is also mani- 
fested in the low amplitude of their two-point correlation 
function relative to that of the overall halo population; see 
Figs. ||, ||). This result can be interpreted as correlation 
between the time since the last major merger and present- 
day environment of the halo (i.e., halos which underwent 
an early major merger tend to be located in clusters and 
groups). The significance of the results was estimated by 
a non-parametric Monte Carlo test which showed that 
the segregation and anti-segregation have significance of 
> 95% in the distance range 0.8/i _1 Mpc to 2.8/i~ 1 Mpc . 
Similarly, the probability of finding stripped halos in pairs 
of separations < 5/i~ 1 Mpc is twice higher than the cor- 
responding probability for the overall halo sample. Halos 
which experienced early major mergers and/or mass loss 
due to tidal stripping are likely to host early type galax- 
ies. In this case, the above mark correlation results for DM 
halos indicate that morphological segregation of galaxies 
may be due to the specifics of the mass evolution histories 
and environment of their parent halos. 

The analysis presented in this paper showed that 
MCFs provide powerful, yet algorithmically simple, quan- 
titative measures of segregation in halo spatial distribu- 
tion with respect to their properties (e.g., maximum cir- 
cular velocity) and merger history (e.g., time since the last 
major merger). The mark correlation functions allow us to 
quantify the degree of segregation as a function of scale 
and can be used to quantify the differences in the spa- 
tial distributions of various galaxy samples (similarly to 
the usual two-point correlation function) and, at the same 
time, to study the interplay between the spatial cluster- 
ing and the distribution of galaxy properties (marks) . In 
this respect, the MCFs are a natural extension of the spa- 
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